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O ■ Abstract 

a^ 

' We present measurements of triple gauge boson coupling parameters using data recorded by 
1^ the OPAL detector at LEP2 at a centre-of-mass energy of 172 GeV. A total of 120W-pair 



P-i! candidates has been selected in the qqqq, q^iVi and iViJ'ui' decay channels, for an integrated 
^ I luminosity of 10.4 pb~^. We use these data to determine several different anomalous coupling 
parameters using the measured cross-section and the distributions of kinematic variables. We 
- measure aB<^=0.35ll:g? ± 0.38, avF^=0.00l[5:i ± 0.11, avK=0.18tH? ± 0.23, Ac/^=-0.03l[5:3? ± 
^ : 0.14, AK\^^^^'>=Omtofi ± 0.20, and AK=0.03l[|i^ ± 0-21. Combining the aw^ resuh with our 
5^ , previous result obtained from the 161 GeV data sample we measure q;vi/</)=— O.OStgls ±0.10. 
All of these measurements are consistent with the Standard Model. 
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1 Introduction 



During the initial year of operation of LEP2, centre-of-mass energies have been attained which 
allow the production of W^^'W^ boson pairs for the first time in e"'"e~ collisions. A total 
integrated luminosity of 9.9 pb^^ was recorded at a centre-of-mass energy of 161 GeV and 
10.4 pb"^ at 172 GeV. 

The W~*'W~ production process involves the triple gauge boson vertices between the W"'"W~ 
and the T? or photon. The measurement of these triple gauge boson couplings (TGCs) and the 
search for possible anomalous values is one of the principal physics goals at LEP2. Previous 
direct measurements have been presented by CDF |jl[ and DO p[, and by the LEP experiments 
using the 161 GeV ^ and 172 GeV p data. In this paper we present an analysis of the 172 
GeV data sample to make measurements of several anomalous coupling parameters. We also 
combine these results with the measurement obtained from our 161 GeV data sample 0. 

The observable effects of possible deviations from the Standard Model have been studied ex- 
tensively 0. Anomalous TGCs can affect both the total production cross-section and the shape 
of the differential cross-section as a function of the W~ production angle. Additionally, the 
relative contributions of each helicity state of the W bosons are changed, which in turn affects 
the distributions of their decay products. In this analysis we make use of both production rate 
and angular distributions measured using events where the W boson pair decays to the qq^F^ 
channel (hereafter qqfl^^ also refers to the charge conjugate state). We also use the measured 
production rate for events where the W boson pair decays to the qqqq and lV(l v^i channels. 



Anomalous coupling parameters 

The most general Lorentz invariant Lagrangian |]^, |^, ^, Q which describes the triple gauge 
boson interaction has fourteen independent terms, seven describing the WW7 vertex and seven 
describing the WWZ vertex. This parameter space is very large, and it is not currently possible 
to measure all fourteen couplings independently. Assuming electromagnetic gauge invariance 
and C and P conservation the number of parameters reduces to five, which can be taken as 
g\, Kz, K^, Az and 01- the Standard Model gl = 1^^ = = 1 and A^ = A^ = 0. 

Different sets of parameters have also been proposed which are motivated by SU(2)xU(l) 
gauge invariance and constraints arising from precise measurements at LEPl. One such set 
proposes three independent linear combinations of the couplings |]^, |^ which are not tightly 
constrained [l^, [ll|] by existing LEPl data. These are: 

aB,p = Ak^ - A^if cos^ Ou, (1) 
aw<f, = (2) 
aw = A^. (3) 

with the constraints that Ak^ = — A/t^tan^ 6^ + Agl and A^ = A^. The A indicates the 
deviation of the respective quantity from its Standard Model value and 6yj is the weak mixing 
angle. Each of the a parameters has the value zero in the Standard Model. A similar approach. 



HISZ [|TT|, reduces this set to two parameters with the extra constraint aB(j, = Oiw<i>- This extra 
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constraint is equivalent to Agl = Ak^/(2cos^ 6^), and Ak,^ and A-y = aw are normally used as 
the variable parameters. 

In this analysis, we present measurements of the following parameters: 

• astj), Q-wtf) and aw, 

• the two and three dimensional correlations of aBip, aw<p and aw, 

• Agl and the two dimensional correlation of Agl and aw, 

• Ak,^ assuming the HISZ constraints, and the two dimensional correlation of A/t^ and aw, 

• Ak, assuming Ak^ = Ak^, and the two dimensional correlation of An and aw- 

In each case, all parameters not mentioned explicitly are set to their Standard Model values. 
The measurements of Ak^ and Ak are included to facilitate comparisons with results from the 
CDF and DO collaborations. 



2 Data selection and Reconstruction 



The data were recorded at an average centre-of-mass energy |jl2[ of ^/s = 172.12 ± 0.06 GeV 



using the OPAL detector, which is fully described elsewhere |T^, A total integrated lumi- 
nosity of 10.36 ± 0.06 pb^^ was recorded at this energy. Events are selected corresponding to 
different W pair decay combinations. These are: (i) qqqq events where both W bosons decay to 
a quark-antiquark final state, (ii) qjqiVg events where one W decays to a quark-antiquark final 
state, and the other W decays to an electron, muon or tau, plus a neutrino and (iii) iu^i v^i 
events where both W bosons decay leptonically. 

The selection procedures which have been developed for these channels are based upon a multi- 



variate likelihood method. These procedures are described fully in |15[ and the results presented 



there are summarised in table |^. In this analysis we use all of the efficiencies, background and 



systematic error evaluations exactly as presented in [|I^ to determine the constraint upon 



anomalous couplings arising from the observed event rate. 

We also constrain anomalous couplings using the production and decay angular distributions in 
the Gfitui channels. We choose to use only these channels as they are the most straightforward 
to reconstruct and there is no ambiguity in assigning decay fermion pairs to each W, nor in 
determining the charge of each W through the decay lepton charge. For this part of the analysis 
we make additional requirements to those described in |T5| in order to ensure that the charged 
lepton is well identified and to further reduce background contamination. In the following we 
describe the reconstruction of qq^F^ events used for the angular distribution analysis and the 
extra selection requirements and the corresponding reevaluation of the backgrounds. 

All Monte Carlo samples which are used are generated assuming the current central value of 
the world average W boson mass of Mw = 80.33 ± 0.15 GeV |]16|. The total W"'"W~ cross 
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section is assumed to be 12.4 pb as obtained from the GENTLE [11] program. All Monte Carlo 
samples used are passed through the full OPAL simulation program [|l^] and then subjected to 
the same reconstruction procedures as applied to the data. 



2.1 qqeL/g and qc^fiu^ events 

The qqeI7g and oslfiT'^ events are characterised by two well-separated hadronic jets, a high mo- 
mentum charged lepton and missing momentum from the unobserved neutrino. The selections 



described in [|15| result in a single track being identified as the most likely lepton candidate. 
The electron momentum is constructed from the direction measured by the tracking detectors 
and the energy measured in the electromagnetic calorimeters. In the case of muons, the mo- 
mentum measured using the tracking detectors is used. The remaining tracks and calorimeter 



clusters in the event are grouped into two jets using the algorithm |T9[. The total energy 
and momentum of each of the jets are calculated using both tracks and electromagnetic clusters 
using the method of pO[. The main residual background is from misidentified qqrz7r events. 



We reduce this further by demanding M^;^ > 50 GeV, where M^i, is the invariant mass con- 
structed from the lepton four-momentum and the missing three-momentum assuming a zero 
mass. This extra requirement results in efficiencies of 85% and 87% for the electron and muon 



channels respectively, which is a loss of 0.5% in each channel with respect to reference |jT^. The 
expected number of signal events and the number of events selected from the data are shown 
in table |l[ The values in this table correspond to an integrated luminosity of 10.36 ± 0.06 pb^^ 
at 172.12 ± 0.06 GeV. Table ^ shows the signal and background cross sections used in the 
analysis of the angular distributions. The expected number of signal events is evaluated using 
PYTHIA [pTI and is normalised such that the total cross section agrees with GENTLE. 



2.2 qqTL^r events 

The signature of qqrz/^ events is less well defined than that of the electron or muon channels 
due to the presence of one or more additional neutrinos, and the absence of a high momentum 



lepton. The primary selection for this channel described in results in the identification of 
the most likely one or three track tau decay candidate classified as an electron, muon, one-prong 
hadronic or three-prong hadronic decay. The remaining tracks and calorimeter clusters in the 
event are grouped into two jets as described above. 

The ensuing analysis of this channel relies upon the correct tau decay products being identified 
in order that they give a reasonable approximation to the original tau fiight direction, and that 
the direction of the hadronic W system can be reliably measured. We therefore make additional 
requirements. In the electron and muon decay channels we require the lepton momentum to be 
greater than 5 GeV. We demand that the track is separated from other activity in the detector 
by requiring that the angle between the lepton direction and both jet axes be greater than 20°. 
In the one-prong and three-prong hadronic decay channels we require the total momentum of 
the decay products to be greater than 11 GeV and the angle between the direction of the vector 
sum of the decay products and either jet axis to be greater than 20°. In addition we require 
the two hadronic jets to be back-to-back within 65°. These requirements result in a correct 



6 



tau decay product identification of 97%, 96%, 90% and 89% in tlie electron, muon, one-prong 
liadronic and tfiree-prong liadronic decay cliannels respectively, where the tau decay products 
are said to be correct if they lie within 10° of the original tau flight direction and have the 
correct charge. The overall selection efficiency drops from 61% to 53% with respect to |jl5|. 
The expected signal and the total number of events selected from the data are shown in tables ^ 
and |. 



2.3 Backgrounds 



The estimated background contaminations from different sources are given in table |^. The 
main background process, Z°/7 qq, is simulated using PYTHIA. Most of the processes lead- 
ing to four fermions in the final state are evaluated using grc4f and EXCALIBUR p3 



Backgrounds from two-photon processes are evaluated using PYTHIA. 

The background labelled in table ^ as '4-fermion' represents the contribution due to all four 
fermion diagrams excluding the three main W+W^ production diagrams (known as "CC03" 0). 
This contribution includes all diagrams which lead to the same final state as the signal, and 
which may therefore interfere to give a negative contribution, as well as those leading to dif- 
ferent final states but which can be misidentified as signal. The 4-fermion contamination is 
obtained using grc4f, by subtracting the accepted cross section for events generated using only 
the CC03 diagrams from that found for events generated using all diagrams. The results were 
found to be in good agreement with a similar study using EXCALIBUR. 



Channel 


Expected Signal 


qqiuf, Bkg. 


Other Bkg. 


Total 


Observed 


From |15| 


[see note in caption 




W+W-- 


' qqeI7e 


16.2 ±0.4 


0.6 ±0.0 


1.2 ±0.4 


18.0 ±0.6 


19 


W+W-- 


> qq/iI7^ 


16.5 ±0.4 


0.9 ±0.0 


0.6 ±0.1 


17.9 ±0.4 


16 


W+W-- 


> qqrur 


11.6 ±0.4 


1.9 ±0.1 


2.7 ±0.7 


16.2 ±0.8 


20 


W+W-- 


y qqqq 


46.8 ± 1.2 




14.3 ±2.9 


61.0 ±3.1 


57 


W+W-- 




10.7 ±0.3 




0.8 


±0.2 


11.5 ±0.4 


8 


After additional requirements for angular distribution analysis. 


W+W-- 


y qqeFe 


16.0 ±0.4 


0.5 ±0.0 


0.6 ±0.3 


17.1 ±0.5 


17 


W+W-- 




16.4 ±0.4 


0.6 ±0.0 


0.2 ±0.3 


17.2 ±0.5 


15 




y qqrI7^ 


9.8 ±0.4 


0.8 ±0.1 


1.6 ±0.3 


12.2 ±0.5 


16 



Table 1: Observed number of candidate events in each W"''W~ decay channel, together with 
expected numbers of signal and background events predicted by the Standard Model. The 
predicted numbers for signal include systematic uncertainties from the efficiency, luminosity, 
beam energy, W"'"W~ cross-section and Mw, while the background estimates include selection 
and luminosity uncertainties. (Note: In this table, contamination in the o^ii'e channels from 
other qq£X/£ channels is listed separately, whereas in |jl5| it is included in the signal.) 
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Channel 


Source 


qqei/e 


qq/^'^M 




Expected signal (fb) 


1543 ± 16 


1579 ± 13 


962 ± 19 


Background (fb) 








W+W~^ qqeFe 




5± 1 


37 ±3 




0.8 ±0.5 




44 ±4 


W+W~^ qqrz7^ 


49 ±4 


59 ±4 




W+W"^ qqqq 


0±0 


0.6 ±0.4 


8±2 




0.8 ±0.5 


0.3 ±0.3 


0.6 ±0.4 


zVt qq 


22 ±6 


9±4 


43 ±8 


4-fermion (non-qqee) 


-1 ± 25 


15 ±25 


59 ±20 


4-fermion (qqee) 


20± 11 


0±0 


44 ±20 


Two-photon 


22± 13 


0±0 


0±0 


Total (fb) 


113 ±31 


88 ±25 


236 ± 30 


Total (%) 


7±2 


6±2 


25 ±3 



Table 2: Signal and background cross sections for the qqez/g, qqAii^^ and qqrI7^ samples used in 
the analysis of angular distributions. The errors on the expected signal include both statistical 
and systematic errors. The errors on the background estimations are statistical only. 

3 Kinematic variables for the qjqiUi event sample 

For each qqiVi event we measure three angles ||, |2^ : 



1) cos6'-w, the production angle of the W with respect to the e beam direction, 

2) cos 6^, the polar angle of the charged lepton with respect to the parent W flight direction 
measured in the W rest frame, 

3) 01, the azimuthal angle of the charged lepton around the axis given by the parent W 
flight direction, measured in the W rest frame. 



In the case of the qqeFg and qq/iz^^i events we perform a three constraint kinematic flt, which 
demands energy and momentum conservation assuming a zero mass for the missing neutrino. 
The flt minimises a sum which also includes a contribution from the difference between the 
nominal world average W mass of 80.33 GeV and the fltted masses of both the hadronic and 
leptonic systems. We demand that the kinematic flt converges with a probability of > 10~^. 
For the small number of approximately 5% of events which fail at this point we revert to 
using the results of a flt without the W mass constraints. We then obtain cos^w by adding 



^The definition of the axes is such that z is along the parent W flight direction and y is in the direction 
X W where e~ is the electron beam direction and W is the parent W flight direction. The axes are defined 
in the W rest frame. 

^ The W mass is treated as Gaussian in the kinematic fit. However, in order to simulate the expected 
Breit-Wigner form of the W mass spectrum, the variance of the Gaussian is updated at each iteration of the 
kinematic fit in such a way that the probabilities of observing the current fitted W mass are equal whether 
calculated using the Gaussian distribution or using a simple Breit-Wigner. 
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together the kinematically fitted four-momenta of the two jets. The charges of the W bosons 
are determined by the sign of the charged lepton. The decay angles are obtained from the fitted 
charged lepton four-momentum after boosting back to the parent W rest frame using the fitted 
W four-momentum to determine the boost. 

In the case of the c\(]tT't events we obtain cos 6*-^ by adding together the measured four-momenta 
of the two reconstructed jets, and the charges of the W bosons are determined by the sign of the 
sum of the charges of the tracks of the tau decay products. In order to reconstruct the decay 
angles, the flight direction of the tau is approximated by the direction of its charged decay 
products. The four unknown quantities can then be calculated using energy and momentum 
conservation. These are the energy of the tau and the three-momentum of the tau neutrino 
originating directly from the W decay. The decay angles are then obtained as for the other 
qqiue events. 

In figure |I] we show the distributions of all three angles obtained from the combined qqe!7e, 
qq/iz/^ and qqTVr event sample, and the expected distributions for aw<t> = ±2 and 0. The 
shapes of the the distributions are obtained from fully simulated Monte Carlo event samples. 
The distributions are normalised to number of events observed in the data. 



4 Anomalous coupling analysis 

In this section we present the analysis of the W^W~ event sample in order to place limits upon 
anomalous coupling parameters. The analysis uses the event sample in two distinct ways. The 
first part uses only the qqiVi events to constrain TGCs through the predicted variation of the 
differential distribution of the production and decay angles. We do not distinguish between the 
qqeFe, qq/il^^ and qqrz/^ events. The second part uses all channels to constrain TGCs through 
the variation of the expected cross-section as a function of TGC parameter. A log likelihood, 
logL, curve is calculated for each part of the analysis. These are independent and are added 
together to obtain a combined log L curve from which the results are obtained. In the following 
the symbol a stands for a generic anomalous coupling parameter. 

A binned likelihood method is used to analyse the three dimensional differential cross section. 
We divide cos6'w into 20 bins in the range [—1,1], cos 6*1 into 10 bins in the range [—1,1] and 
(j)} into 5 bins in the range [— 7r,7r]. The logL distribution is obtained in several steps. 

In the first step we parameterise in each bin the expected cross-section due to W^W^ pro- 
duction, before detector and acceptance effects, as a function of a. We denote this as af^'^^a), 
where i refers to each bin in the three dimensional space. It is obtained using large samples of 
EXCALIBUR Monte Carlo events including the effects of the W width, Fw, and initial state 
radiation, ISR, but without detector simulation. In each bin, af^'^la) is parametrised using 
the fact that a is a quadratic function of a (as a enters linearly in the Lagrangian). Several 
samples are generated for different values of a and a quadratic parameterisation obtained to 
predict a{a) for all other values of a. 

In the second step we calculate a correction matrix to include the effects of acceptance, res- 
olution and feedthrough from other W"'"W~ decay channels. We use fully simulated Monte 
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Carlo events generated with a = to obtain factors, c^j, which allow for events generated in 
true bin i being reconstructed in bin k. All W+W" decays which are reconstructed as qq^r'^ 
events in bin k are counted in the correction factors Cki- The c^i therefore include the effect of 
feedthrough from other W"'"W~ decay channels which is always a fixed fraction of the signal 
which is independent of the value of a. For a given bin, corresponding to a limited phase 
space region, the Cki factors are assumed to be independent of a. Possible biases due to this 
assumption are discussed in section ^. 

Combining these terms results in the expected observed cross-section for each bin k, due to all 
W+W~ channels: 



In the next step the cross-section for the Z^/'j — > qq and two photon backgrounds, cr^^^, is esti- 
mated using PYTHIA and is added to give af^{a). The small background due to four fermion 
diagrams given in table |^ is neglected. The effect of this neglect is discussed in section ^. The 
quantity af^{a) now contains all of the information on the shape of the expected distribution 
as a function of a. We multiply this by a scale factor to give the prediction for the number of 
events in each bin, n^, where the scale factor is chosen such that the predicted total number of 
events is equal to the actual number of events observed in the data. This is to ensure that we 
do not incorporate any information from the overall production rate in this part of the analysis. 
The probability for observing the number of events seen in each bin for an expectation of Uk 
is calculated using Poisson statistics. The statistical fluctuations in Cki and a^^^ are taken into 



account using the method of reference |25]. The negative logL distributions for each of the a 
parameters studied are shown in figure |^. 

The information given by the observed event rates is included by calculating the likelihood for 
the mean number of events expected as a function of a to have resulted in the observed number 
of events given in the upper part of table |T]. This is calculated separately for each channel. 
The variation of the mean number of signal events expected as a function of a is determined 
using fully simulated EXCALIBUR events subjected to the same selection requirements as 
the data. The number of selected events is parameterised using the quadratic dependence 
mentioned earlier. The overall normalisation is adjusted to agree with the GENTLE prediction 
at the Standard Model point. The probability of finding the observed number of events is 
calculated assuming a Poisson distribution for the signal and background in each channel. 
The probabilities obtained for each channel are multiplied together and the negative logL 
distribution obtained. These are shown in figure |^ for each parameter considered. 

The curves due to the total and differential cross-section information are consistent with each 
other for each of the a parameters. Both curves provide a useful constraint. These are inde- 
pendent measurements and are added together to give the overall log L distributions shown as 
the dash-dotted lines in figure 

We also calculate the two dimensional likelihood distributions resulting when two anomalous 
couplings are simultaneously allowed to vary from the Standard Model. This is incorporated 
into both parts of the analysis described above by parameterising the predicted variation of 
signal events using a quadratic function of both TGCs. The resulting 95% probability contours 
for different pairs of couplings are shown in figure |^. We also perform the analysis allowing 
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all three of the a parameters to vary simultaneously. The resulting set of two dimensional 
correlation contours are also shown in figure ^ a), b) and c), where in each case the third 
parameter is varied in order to re- minimise the logL at each point in the two dimensional 
space. 



Alternative analyses 



We have developed two other independent analyses, which are largely complementary. The 
first method uses a simple unbinned maximum likelihood approach which does not include 
the effects of Fw, ISR, acceptance, resolution or background in the analysis of the kinematic 
variables. These effects are expected to be small in comparison to the present level of statistical 
precision, and so the method provides a robust systematic comparison. The second method 
uses an optimal observables method |2^, which includes all detector effects. The log L curves 
produced by all three methods agree very well and give us confidence in the stability of the 
results. 



5 Systematic error studies 



The Monte Carlo simulation of the measured quantities depends mainly upon the simulation of 
the jets from the W hadronic decay. Jet reconstruction has been studied and tuned extensively 
at LEPl, showing good agreement between distributions measured from data and Monte Carlo 
samples. We therefore expect this to be adequate for the small number of events in this sample. 
Studies of back-to-back jet pairs using the LEPl data yield the following possible differences 
between data and the Monte Carlo simulation of jets; 10% for energy resolution, 0.5% for energy 
scale and 10% for resolution in cos^ and 0. As a direct test of cos6'w reconstruction we have 
used radiative Z° qq events taken from the = 91, 161 and 172 GeV data. By selecting 
events containing observed radiated photons with energies up to 20 GeV we obtain a sample 
of jet pairs exhibiting a similar acoUinearity distribution to W decays. Assuming that no other 
photons have been radiated in the event the true direction of the Z° qq system is opposite 
to that of the photon and the difference between the value of cos Q measured from the photon 
and that measured from the hadronic system is therefore strongly related to the resolution of 
cos^w The results obtained from both data and Monte Carlo events agree well showing no 
significant differences in shape or width of the distributions. We obtain a conservative upper 
limit upon the possible relative shift on cos^w between data and Monte Carlo of 0.01. All of 
these variations are then used to vary the Monte Carlo jet reconstruction in the analysis and 
the resulting changes caused to the TGCs were added in quadrature and taken as a systematic 
error. These are shown in table ^ labelled 

We have evaluated the systematic errors due to several other factors. 



b) The effect of limited Monte Carlo statistics. 

c) The uncertainty in the world average measured W mass and the LEP beam energy. 

d) The uncertainties of the overall normalisation (selection efficiency and luminosity) given 



m 



15| used in the calculation of the logL curve derived from the event rate. 
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e) We have varied the Monte Carlo generator to use both EXCALIBUR and KORALW 



to generate the Standard Model sample. We have varied the fragmentation model by 



using HERWIG ^ 



m 



f) We have varied the backgrounds within the estimated systematic errors given in 
the calculation of the log L curve derived from the observed event rate. We have doubled 
and removed the TP /'^ and two-photon background additions used in the calculation of 
the log L curve derived from the differential distributions. 

g) The analysis was carried out using Monte Carlo samples generated at a nominal energy 
of 171 GeV. We have estimated the possible bias this introduces by using Monte Carlo 
samples generated at 172 GeV as test data. 



The errors due to these sources are listed in table 0. 

This analysis relies upon EXCALIBUR correctly describing the variation of cross section, an- 
gular distribution and four-fermion background as functions of anomalous TGCs. We have 
investigated the limits of this assumption by comparing these quantities to those predicted 



by several other programs (GENTLE, KORALW and ERATO Q). We find that the total 
cross section predictions agree to within approximately 2.5% across a wide range of anomalous 
TGC values and we add this variation as a further systematic error to d). We find that both 
the angular distributions and the contributions of diagrams in addition to the three W"'"W~ 
diagrams are also compatible between the programs at the level of a few percent and we assign 
no extra systematic error. 

Finally we consider possible biases of the method. This method allows for all effects of Fw, 
ISR, acceptance and detector resolution in the calculation of both log L curves and is therefore 
in principle free from biases due to these effects. However biases may be introduced by the 
coarse binning, by the approximation in the correction used to include detector effects, or 
by the neglect of the small four fermion background. To quantify these effects we perform 
the analysis using samples of approximately 20,000 fully simulated four fermion Monte Carlo 
events produced by both grc4f and EXCALIBUR as input. These samples were generated with 
different TGC values. In all cases the TGC value reconstructed by the analysis is found to be 
consistent with the generated value and no significant difference was observed between the grc4f 
and EXCALIBUR samples. The greater of the statistical precision of the test or the measured 
bias was taken as a conservative systematic error and is shown in row h) of table |^. 

All contributions to the systematic errors are summarised in table |^. To include the systematic 
error in the likelihood curves, we apply a simple procedure which would be exact if all errors 
were Gaussian. We calculate a function of a, /S.\og L{a)syst, to add to the logL curves. This is 
given by 

A\ogL{a)syst = {oi- aoYliVstat + Kyst)"^ - {VstatY\a - ao) 

where a represents a vector of one or more TGC parameters, ao is the TGC value at the 
minimum of the logL curve, Vgyst is the systematic error matrix constructed assuming the 
errors are Gaussian, and Vstat is approximated (for this purpose only) by assuming a parabolic 
logL distribution around Oq- The full logL curves for each parameter are obtained by adding 
A log L{a)syst to the log L curve due to statistical errors only, and these are shown as solid lines 
in figures and ^ 
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Source 




Error on 


Darameter 








OiB(f> 


aw<j> 


aw 


Agl 




Ak 


a) 


Jet reconstruction 


0.21 


0.05 


0.03 


0.05 


0.12 


0.12 


b) 


MC statistics 


0.06 


0.04 


0.05 


0.05 


0.06 


0.07 


c) 


W mass/LEP energy 


0.04 


0.01 


0.01 


0.00 


0.01 


0.01 


d) 


Overall normalisation 


0.07 


0.01 


0.02 


0.01 


0.02 


0.02 


e) 


MC generator 


0.11 


0.02 


0.15 


0.02 


0.04 


0.08 


f) 


Background 


0.22 


0.06 


0.12 


0.06 


0.11 


0.09 


g) 


171/172 GeV diff. 


0.12 


0.05 


0.10 


0.08 


0.07 


0.07 


h) 


Bias tests 


0.11 


0.04 


0.05 


0.05 


0.06 


0.06 




Total 


0.38 


0.11 


0.23 


0.14 


0.20 


0.21 



Table 3: Contributions to the systematic errors in the determination of TGC parameters. 

6 Results and Conclusion 

We express the results as measurements with the one standard deviation limits given by the 
values of a where the negative logL curve rises by 0.5 from the minimum. The results are 
given in table ^ where the statistical and systematic errors are shown separately. The expected 
statistical errors for a sample of this size were evaluated by repeating the analysis using Monte 
Carlo events divided into sub-samples equivalent to the actual data sample. The RMS spread 
and the average error agreed within 10%, and the latter is shown in the table. 

We also express the results as 95% confidence level (C.L.) limits given by the values of a where 
the change in likelihood is 1.92. The ranges obtained are shown in table ^. 

In our previous publication we have presented a measurement of the aw,f, parameter using 
the 161 GeV data. By including the logL obtained from that measurement in this analysis 
we obtain the combined result aw,f)=—^-^^^o'M> 0.10. The corresponding 95% C.L. range is 
-0.59 < aw<t> < 0.52. 

The results which we have obtained for the parameters listed as c), e) and f) in the results 
table can be compared with those reported recently by CDF |Q and DO which are a factor 
of between two and three more precise. The measurements agree in all cases. 

In conclusion, we have presented results using W"'"W~ events reconstructed from 10.4 pb~^ of 
data recorded at LEP2 at a centre-of-mass energy of 172 GeV. We have used the observed 
W"'"W~ event rate and the distribution of kinematic variables to place limits upon possible 
anomalous triple-gauge-boson couplings. We measure a_B0=O.35li;o7 ± 0.38, avy0=O.OOlo!28 i 
0.11, avF=0.18l|]:^?±0.23, A^^=-0.03i|]:^?±0.14, A4^^^^)=0.03l[5if ±0.20, and A/t=0.03l[;i^± 
0.21. These results agree with the Standard Model expectation of zero for each anomalous cou- 
pling. The negative logL distributions for these parameters are shown in figure ^ The two 
dimensional correlation contours are shown in figure ^. 
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Parameter 


Measurement 


Expected 
stat. err. 


95% C.L. 


range 


OiB(t> 


0.35li;o7 ± 0.38 


±1.1 


r -1 

[-1.6, 


2.7] 


aw<f> 


O.OO^ais ± 0.11 


±0.29 


[-0.55, 


0.64] 


aw 


0.18lo:47 ± 0.23 


±0.48 


[-0.78, 


1.19] 




-o.o3l[]:;^?±o.i4 


±0.39 


[-0.75, 


0.77] 




o.o3lg:i? ± 0.20 


±0.53 


[-0.98, 


1.24] 




0.031^6 ± 0.21 


±0.50 


[-0.90, 


1.12] 



Table 4: Measurements of anomalous coupling parameters using the OPAL 172 GeV data 
sample. The measured values and one standard deviation errors are shown in column 2, where 
the statistical (asymmetric) and systematic (symmetric) errors are shown separately. The 
expected error obtained from multiple Monte Carlo samples is shown in the third column. The 
95% confidence level ranges including systematic errors are shown in the fourth column. 
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Figure 1: The distributions of the kinematic variables cos6'w, cos 9}, (J)} obtained from the 
qqii^i events. The hatched histogram shows the non-qq^I/^ background. These are compared 
with the distribution expected in the Standard Model using fully simulated Monte Carlo events. 
The predicted distributions for aw<t> = +2(— 2) are also shown as dotted (dashed) lines. Note: 
in the case of W"*" Iv events the value of ^| is shifted by tt in order to overlay W"*" and W~ 
distributions in the same plot. 
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OPAL 




Figure 2: The negative log likelihood curves (statistical contribution only) obtained from 
the total event rate (dotted) and the shape of the differential distributions (dashed). The 
dot-dashed curve is obtained by adding these together. The solid curve shows the final result 
obtained when the systematic error is added. The A prefix indicates that the likelihood is 
shown relative to its minimum value. The figures correspond to a) ccs^, b) aw(i>, c) aw, d) Agf, 
e) Ak-y assuming the HISZ constraints and f) Ak assuming Ak^ — Akz. 
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OPAL 




Figure 3: The 95% CL. two dimensional correlation contours for different TGC parameters. 
The effect of statistical and systematic errors is included. The star indicates the minimum 
point. The first three plots show all pairs of the cts^, aw4>, and aw parameters. Figure d) 
shows Agl-vs-aw, figure e) shows AK^-vs-aw assuming the HISZ constraints and figure f ) shows 
An-vs-aw assuming An^ = Ak^- In figures a), b) and c) the dashed line shows the results 
obtained by varying all three a parameters (see text), where both statistical and systematic 
errors are included. 
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